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The neutron-rich mass region around A ~ 100 presents challenges for modeling the astrophysical 
r-process because of rapid shape transitions. We report on mass measurements using the TITAN 
Penning trap at TRIUMF-ISAC to attain more reliable theoretical predictions of r-process nucle- 
osynthesis paths in this region. A new approach using highly charged {q — 15+) ions has been 
applied which considerably saves measurement time and preserves accuracy. New mass measure- 



ments of neutron-rich s^.s'i'.ssj^^ ^^^^ 



Sr have uncertainties of less than 4 keV and show 



deviations of up to llcr to previous measurements. An analysis using a parameterized r-process 
model is performed and shows that mass uncertainties for the A = 90 abundance region are elimi- 
nated. 

PACS numbers: 21.10.Dr, 26.30.Hj, 27.60. +j, 82.80.Qx 



I. INTRODUCTION 

The atomic mass, and from it the derived nuclear 
binding energy, is a key property of the nuclear many- 
body system. The binding energy reflects all interactions 
among the constituents and provides important informa- 
tion for a broad variety of studies in nuclear and atomic 
physics [l|, [2| including tests of the standard model of 
particle physics jll, nuclear structure studies and tests 
of theory and models (e.g. f4]), and nuclear astrophysics 
0, including nucleosynthesis pathways. In addition, the 
conserved vector current (CVC) hypothesis and the uni- 
tarity of the Cabbibo-Kobayashi-Maskavifa (CKM) quark 
mixing matrix ([si] and references therein) are probed by 
determining Vud from super-allowed /3 emitters. One of 
the key inputs is the Q-value obtained from the mass 
difference between the mother and daughter of the de- 
cay. Nuclear structure variations from nucleus to nucleus 
can be sensitively probed by differences in nuclear masses 
and thus provide insight to sudden or gradual changes in 
nuclear structure. (See discussion in Q, recently con- 
firmed by 7-ray spectroscopy by Q-) Here we report on 
new measurements of masses of neutron-rich nuclei in the 
A w 100 region. 

The synthesis of about half of the heavy elements be- 
yond germanium (Z = 32) proceeds in nature via the 
rapid neutron-capture process (r-process) Q. In the 
most common r-process models the neutron capture reac- 
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tions proceed until an equilibrium between neutron cap- 
ture and photo disintegration, (n,7) < — > (7,n), is estab- 
lished, driving the r-process path to nuclei with neutron 
separation energies of approximately 2 to 3 MeV Q. 

Testing r-process models against abundance observa- 
tions requires reliable nuclear physics that translates a 
set of evolving astrophysical conditions into a character- 
istic nuclear abundance pattern. Nuclear input data of 
importance to r-process models are nuclear masses, /3- 
decay half-lives and branching ratios for neutron emis- 
sion. Fission rates and fission fragment distributions, 
neutrino interaction rates, and charged particle fusion 
rates also play a role Q. 

In models characterized by an (n,7) i — >■ (7,n) equilib- 
rium, the reaction path for a given set of astrophysical 
conditions is governed by nuclear masses. However, many 
of the nuclei in the r-process are inaccessible experimen- 
tally. Hence, nuclear astrophysics calculations rely heav- 
ily on theoretical mass predictions, which are performed 
by models with parameters that are adjusted to known 
masses [Ij. Experimental mass data on neutron-rich nu- 
clei are therefore not only needed as direct input into 
r-process model calculations, but also to improve mass 
models and reduce the need for extrapolations. Fig. [1] 
shows the reaction flows of the parameterized, fully dy- 
namical r-process model that is used in this work, fol- 
lowing Freiburghaus et al. ;10] and Hosmcr et al. [T3j. 
The model is based on an adiabatic expansion as it might 
be encountered in high entropy neutrino driven winds in 
core collapse supernovae. The path is characterized by 
a complex network of charged particle and neutron in- 
duced reactions and their inverse. More details can be 
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FIG. 1. (color online) A section of the nuclide chart. The neutron-rich nuclei studied in this article are indicated by red stars. 
The mass accuracy is displayed as well as the calculated time-integrated net-reaction-flows for a single S = 100 component 
(blue) of the parameterized, high entropy wind inspired, r-process model used in this work. Flows above a relative final 
abundance of 10""' are denoted by light blue lines indicating the complex interplay of charged-particle and neutron-induced 
reactions in this mass region; dark blue lines mark the outer boundary of the reaction flows for clarity. 



found later in the text. The model passes through the 
mass region covered by our experiment {A « 100), which 
helps in reducing the challenge for theoretical mass mod- 
els. 

Additional motivation for this work stems from the de- 
sire to better understand the nuclear structure, in partic- 
ular shell closure effects at = 50, the subshell closure 
at A^ = 54, and a possible onset of large deformation for 
nuclei with > 60 ^2i|. Several theoretical investiga- 
tions of the nuclear structure in the N ^ 60 region have 
been carried out ([Tsl-fTsI and references therein). Within 
a self-consistent Hartree-Fock-Bogoliubov (HFB) formal- 
ism [l3|, the nearly-spherical shape for A^ < 59 Sr, Zr, 
and Mo isotopes is predicted to develop into a strongly 
deformed prolate shape above A^ > 60. In the so- 
called DlS-Gogny energy-density-functional framework 
( [l5| and references therein), in which one-quasiparticle 
configurations are employed, a self-consistent mean-field 
approximation indicates nuclear deformation and shape 
co-existence. These calculations predict a sharp oblate- 
to-prolate transition in the Rb, Sr, and Zr isotopes and 
triaxiality in the Mo isotopes. 

Empirical evidence supports rapidly changing behavior 
in A^ 60 nuclei @, 0, [l6l - [20| . For example, nuclear de- 
formation can manifest itself in the two-neutron separa- 
tion energy S2n, the difference in energy in isotopes sepa- 
rated by two neutrons. Typically, S2n decreases smoothly 
with the neutron number N, and a change in slope may 
indicate the onset of deformation or (sub)shell closures. 
This signature has already been observed in the Rb iso- 
topes [2l|, and in the Sr and Zr isotopic chains around 



N fa 60 [3]; however, neutron-rich Kr isotopes do not 
present any nuclear quantum phase transition [6| . In this 
work, we extend the investigation of nuclear deformation 
around A^ 60 to neutron-rich Rb and Sr isotopes via 
Penning-trap mass spectrometry. 

Previous Penning trap mass measurements in this re- 
gion [l^, ^3 differ from the atomic mass evaluation 
(AME03) [il by up to 300 keV and up to llcr. Therefore, 
an independent confirmation was desired. To that end, 
precise and accurate mass measurements on radioactive, 
short-lived isotopes have been performed at TRIUMF's 
Ion Trap for Atomic and Nuclear science setup (TITAN) 
[23I [131 at the radioactive beam facility for Isotope Sep- 
arator and Accelerator (ISAC) [25|. A unique feature of 
TITAN over other Penning traps at radioactive ion beam 
facilities is the ability to charge-breed ions prior to the 
mass measurement. The advantage of high charge states 
is the improved precision, which scales directly with the 
charge state q. Nuclei with short life times can benefit 
from this approach because it enables high precision with 
short observation times in the Penning-trap. This has re- 
cently been demonstrated for ^'*Rb 26] {T1/2 = 65 ms). 
As long as the production yields of the nuclei are high 
enough to counteract additional losses introduced by the 
charge breeding process, the full precision gain of a fac- 
tor of q can be exploited. At lower yields the actual gain 
is reduced but will stay above unity provided that the 
overall charge breeding efficiency is greater than l/q^. 
Alternatively, by boosting the precision the required un- 
certainties can be reached in shorter times which enables 
the mapping of larger parts of the nuclear chart during 
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the limited measurement times at rare beam facilities. 
The experimental time available to measure the seven 
masses in this work was only sa40 hours which was suit- 
able to achieve a precision that is relevant for nuclear 
astrophysics as described later in this publication. A 
conservative approach concerning systematic uncertain- 
ties was taken into account. 



II. EXPERIMENTAL SETUP 

TITAN is a multi-trap setup located at the ISAC [H] 
facility at TRIUMF. The radioactive isotopes are pro- 
duced by bombarding an ISOL-type target with 500 MeV 
protons from the cyclotron. For the first time, a VCx tar- 
get was used at TRIUMF and neutron-rich isotopes were 
produced for this experiment by bombarding it with a 
proton current of 2 ^A [23 • The ion beam was generated 
using a surface-ionization source and extracted at 20 keV, 
transported through the ISAC mass separator with a re- 
solving power m/5m sa 3000 and delivered to the TITAN 
facility. TITAN, shown schematically in Fig. [51 presently 
consists of three traps: a buffer-gas-filled radiofrequency 
quadrupole (RFQ) trap [28] for cooling and bunching, 
an electron-beam ion trap (EBIT) 29] for charge breed- 
ing, and a hyperbolic Penning trap (MPET) for high- 
precision mass measurements on short-lived nuclei with 
a precision of down to 5m /m « 10~® (30j . 

While previous work at TITAN focused on mass mea- 
surements of singly charged ions (SCI), in particular 
light neutron-rich isotopes and isotopes with neutron 
halos [3TI - [36t . the present work used charge-bred ions 
(see Fig. Mass measurements of short lived, highly 
charged ions (HCI) were recently performed at TITAN 
[26l [37I for the first time and here we extend this to 
neutron-rich Rb, Sr isotopes. In this scheme the ions are 
first cooled and bunched in the RFQ, and then sent to 
the EBIT (see Fig. [2]) where they are recaptured. An en- 
ergetic electron beam (up to 70 keV) removes bound elec- 
trons of the initially singly charged ions through impact 
ionization. In this experiment the ions were charge-bred 
for Tbreed ~ 80 ms using an electron beam current and 
energy of 30 mA and 2.5 keV, respectively. Ions were 
then extracted by opening the trap barrier for 800 ns, 
generating a short pulse. This leads to a sharp time sep- 
aration of ions in various charge states or with different 
m/q (see Fig. [3]) and allows for the selection of ions with 
a certain (m/g)-ratio using a Bradbury- Nielsen ion gate 
[28!]. A distribution of charge-bred ®^Rb-ions extracted 
from the EBIT is shown in Fig. [3] In this case a charge- 
breeding time of 197 ms was used. The HCI with the 
charge state taken to be the one with the greatest ratio 
of ion of interest to background are sent to the MPET 
to determine the ion's mass using the time-of-flight ion 
cyclotron resonance (TOF-ICR) method [H, . 

In the MPET ions are trapped radially via a strong 
homogeneous magnetic field B and axially via a weak 
harmonic electric field. The mass of the ion of interest 
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FIG. 2. (color online) The TITAN facility for high-precision 
atomic mass measurements is presently composed of three ion 
traps. The singly charged ion (SCI) beam from ISAC passes 
through an RFQ (radio-frequency quadrupole) where the ions 
are thermalized in a helium buffer gas, bunched, and then sent 
straight to the precision Penning trap for mass measurements 
(MPET), indicated by the dashed red arrow. Alternatively, 
with solid blue arrows the pathway for mass measurement of 
highlycharged ions (HCI) is shown. The SCI from ISAC are 
sent to the EBIT (electron-beam ion trap), charge-bred, and 
transported to MPET. 
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FIG. 3. (color online) Time-of-flight spectra of charge-bred 
ions extracted from the EBIT. A charge-breeding time of 
197 ms, a magnetic field of 3 T, an electron beam current 
and energy of 30 mA and 2.5 keV, respectively, and an ex- 
traction time of 800 ns were applied to obtain this spectra. 
The data represent 500 ion bunches. The red line shows the 
case of injection of ^^Rb^"*" into the EBIT. The black spectrum 
represents background gas in the EBIT. 
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M in the charge state q is determined from the ion's 
cyclotron frequency 



q-B 
27r • M ' 



(1) 



A direct measurement of is done by applying a 
quadrupolar radio-frequency field u-^i which converts the 
radial magnetron motion into the reduced cyclotron mo- 
tion (for details on the application of the technique at 
TITAN see [13] )■ When v-^f is equal to the ions gain 
radial energy, which is converted into longitudinal energy 
when extracted from the trap through the magnetic-field 
gradient and subsequently detected with a microchannel 
plate detector. The resonant ions therefore have a shorter 
time of flight. A frequency range around the expected 
is scanned to obtain a resonance (Fig. |4] and [5|) . 

The statistical uncertainty on a mass measurement is 
given by 



5m 



T ■ 



(2) 



where m is the atomic mass, Trf is the excitation time, 
A'ipn is the number of ions sampled in the measurement 
[ilj , and is a trap specific parameter which is close to 
1 for TITAN. The excitation time Trf is limited by the 
half-life of the radioactive ions. TITAN can boost the 
precision of the mass measurement, which scales linearly 
with q, by increasing the charge state of the ions. Ad- 
ditionally, the increase in resolving power due to charge 
breeding allows the resolution of low-lying nuclear iso- 
mers in Penning traps [i^] . Alternatively, the integrated 
measurement time can be reduced due to the use of HCI 
as shown by this work. For radioactive ions with low 
yields, high charge states allow for a decrease in the to- 
tal number of measured ions while still maintaining high 
precision. In order to assess the benefits of charge breed- 
ing, additional aspects have to be considered. The charge 
breeding process is accompanied by loss mechanisms re- 
sulting from ion transport, capture, storage, and extrac- 
tion. Furthermore, the ions are distributed over several 
charge states, the energy spread is increased, and the 
additional time delay increases decay losses. The prob- 
ability of charge exchange with residual gas during the 
excitation in the Penning trap becomes higher. However, 
for short lived nuclei this excitation time has to be kept 
short. As can be seen in Fig. 2] the quality of the reso- 
nance is not significantly diminished at Trf = 97 ms. As 
demonstrated at SMILETRAP Q with highly charged 
ions of stable nuclei, trapping times of up to several 
seconds are possible without significant charge exchange 
losses. 



III. DATA ANALYSIS 

The main observable in Penning-trap mass spectrome- 
try is the cyclotron frequency i>c (see Eq. [1} . The charge 



state q is derived from a time-of-flight spectrum, but 
the magnetic field needs to be measured. From a fit of 
the theoretical line shape [s^ to the resonance data (see 
Fig. 0] and [5]), the mass can be extracted from Eq. 3] 
if q and B are known. To minimize systematic effects 
(discussed below) and to calibrate the magnetic field, a 
reference measurement of an ion with a well-known mass 
and q is performed before and after the actual measure- 
ment. The primary experimental result is the ratio of the 
cyclotron frequency of the ion of interest to that of the 
reference ion Vc ref , thus the ratio of the masses. 



R = 



Vc q ■ Mycf ' 

The atomic mass m is given by: 



(3) 



m = R ■ (rrirof - Qrcf • rric + Se.ref ) +q-mc- Be, (4) 

where R is the average of all measured frequency ratios. 
Be is the total binding energy of the removed electrons 
(also known as neutralization energy) of the ion of inter- 
est, and the index "ref" refers to the reference ion. 

The value for fc.rof at the time of the measurement of 
Vc is obtained by a linear interpolation of the two refer- 
ence measurements that enclose the measurement of the 
ion of interest. There is a correlation introduced between 
adjacent frequency ratio measurements due to shared ref- 
erences. For the data analysis we take into account a full 
covariance matrix between all the ratios [44] . 

In Tab. U the frequency ratios of 94,97,98j^|-,i5+ g,j^^ 
94,97-99gj,i5+ jgotopes relative to ^^Rb^'^'^ as well as mass 
excesses of the atomic mass are presented. Accuracy lim- 
itations, refiected in systematic uncertainties, such as the 
following are considered: 

• A spatially-uniform magnetic field is required. We 
minimize the effect of instabilities in the mag- 
netic field by using a frequency ratio where sev- 
eral systematics cancel out. As stated in [i^ 
the uncertainty for magnetic field instabilities is 
AR/R <^0.2 ppb/ hour between adjacent reference 
measurements. In this work the time between two 
reference measurements was kept to less than one 
hour. 

• The uncertainty in the mass of the reference ion can 
be considered negligible, since we use ®^Rb which 
is known with a mass uncertainty of 11 eV [45j . 
equivalent to 0.1 ppb. 

• The total electron binding energies were taken from 
[46j . There the total atomic binding energy (bind- 
ing energy of all remaining electrons) has been cal- 
culated using a Dirac-Fock approximation and val- 
ues are tabulated for lithium- to dubnium-like sys- 
tems with Z = 3.. 118. Uncertainties for the total 
electron binding energies of Rb^'^"'"'^^"'" and Sr^^"*" 
are conservatively estimated to be below 20 eV, 
corresponding to 0.2 ppb [47j . 
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TABLE I. Frequency ratios of ^''''■'^'^*Rb^^+ and 94,97-99 gj.i5+ jgQ^opgg relative to *^Rb^^+ as well as mass excesses. Uncertainties 
are expressed in parenthesis. The first error on the frequency ratio represents the statistical error multiplied by the reduced 
chi square of the fit of the line shape, count-class analysis as well as time correlations. The second and third errors represent 
systematics related to A(m/g)-dependent shifts and ambiguities in the choice of time-of- flight range. The fourth error in square 
brackets represents the quadrature sum of all the errors. For the mass excesses the combined uncertainty is shown. In the last 
two columns the mass excess values from AME03 H and JYFLTRAP (Rb [ij, Sr [T!]) are listed. 



Isotope 


7l/2 




A/£titan (keV) 


A/£ameo3 (keV) 


A/i?,]YFLTRAP (keV) 


94j^^l5+ 
97j^^l5+ 

9«Rbi5+ 


2.702 s 
f69.9 ms 
fl4 ms 


0.958672311(22) (16) (3) [27] 
0.989404952(17) (16) (0) [23] 
0.999658513(31) (16) (14) [38] 


-68562.6(2.4) 
-58519.2(2.1) 
-54318.4(3.4) 


-68553.4(8.4) 

-58356.3(30.5) 

-54221.6(50.2) 


-68564(5) 
-58519(6) 


94gj.l5+ 
97s,15+ 
98g^l5+ 
99gj.l5+ 


75.3 s 
429 ms 
653 ms 
269 ms 


0.958559623(10) (16) (0) [19] 
0.989294688(37) (16) (2) [40] 
0.999525849(41) (16) (5) [44] 
1.009776308(42) (16) (3) [45] 


-78845.8(1.7) 
-68581.2(3.6) 
-66424.5(4.0) 
-62506.8(4.1) 


-78840.4(7.2) 
-68788.1(19.2) 
-66645.7(26.3) 
-62185.7(80.0) 


-68587(10) 
-66431(10) 
-62524(7) 



• Although we aim for single-ion injection, the time- 
of-flight spectra of the ions detected on the mi- 
crochannel plate detector after the Penning trap 
show a multitude of (m/q')-states. These are un- 
resolved charge states due to charge exchange with 
residual background gas in the MPET. To minimize 
the possibility of charge exchange a low pressure is 
favored. The pressure in the MPET vacuum section 
was measured to be ss 5 • 10~^^ mbar. The third 
parenthesis in Tab. U represents the uncertainty as- 
sociated with gating on the charge state of interest. 
By varying the time-of-flight range in the analysis 
considering a worst-case scenario, the uncertainty 
is estimated to be a few ppb. 

• If more than one ion is stored in the trap we account 
for shifts due to ion-ion interactions using a count- 
class analysis This uncertainty is convoluted 
with the statistical uncertainty (see first parenthe- 
sis in Tab. |l| and varies from to 3 ppb depending 
on the isotope. 

• Accuracy checks for HCI needed in this work have 
been performed by measuring the cyclotron fre- 
quency of 85Rb"+ vs. 85Rb"+ vs. S5Rbi5+. Dif- 
ferent settings (e.g. Lorentz steerer [i^, extraction 
optics, etc.) as well as different timings (e.g. duty 
cycle, capture timings, etc.) were covered. A shift 
in frequency cannot be excluded due to possible 
trap misalignment, ion-ion interactions due to the 
increase of charge states, and relativistic mass in- 
crease. The mass shift due to the increase of charge 
states scales with the difference in the ratio of m/q 
of ion of interest to reference, A{m/q), the relativis- 
tic effect in this measurement setup with A[q/m). 
The systematic uncertainty for these checks is con- 
servatively estimated to be < 15.5 ppb (absolute) 
and is indicated in Tab. U in the second parenthe- 
sis. Therefore a small A in the ratio of m/q or q/m 
of ion of interest to reference is essential. We used 
®^Rb^'^+ with {m/q) = 6.5 as the reference ion to 



minimize A{m/q) effects. 

• Systematic shifts due to isobaric contamination 
needs to be minimized. The Sr beam was free of 
contaminations whereas in the case of all Rb mass 
measurements, not all contaminating ions could be 
removed using ISAC's mass separator. We used a 
dipole radio-frequency excitation Q preceding the 
quadrupole frequency scan to remove all contami- 
nations. These were identified to be the Sr isobars. 
The application of a dipole (reduced cyclotron) fre- 
quency of the contaminant in the radial plane re- 
sults in the excitation of their cyclotron radii. Dur- 
ing the extraction of the ions the excited ions will 
not clear the aperture. This is a process referred to 
as dipole cleaning. 

IV. RESULTS 

The results of the mass measurements of ^''■^'^■^SRb 
and 94,97-99gj, performed with TITAN are summarized 
in Tab. |T] and discussed in subsections IIV AIIIV Dl The 
absolute uncertainty of all investigated isotopes is below 
4 keV, including the new direct mass measurement of 
^*Rb. In Tab. |T] column Mi^TiTAN hsts the measured 
TITAN mass excess with the different individual uncer- 
tainties, as well as the combined uncertainty. The next 
two columns show results found in the literature for these 
isotopes. The visual comparison between the mass excess 
by TITAN and the previous measurements is shown in 
Fig. [T] A global mass evaluation as outlined in [22*] was 
performed and the evaluated values for the ions of inter- 
est in this work can be seen in Tab. |lT] and later in the 
text. 



The measurement campaign began with ^ Rb 
{Ti/2 =2.702(5) s) which is relatively well known accord- 
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FIG. 4. (color online) ^''Rb^^+ cyclotron resonance taken with 
80 ms charge-breeding time, 20 ms dipole cleaning, and 77 ms 
excitation time in MPET. The solid line is a fit of the theo- 
retical line shape [s^ to the data. 



ing to AME03 [24I, Sm = 8.4 keV. A cyclotron resonance 
of ^^Rb^^"*" is shown in Fig. 21 With our mass excess of 
-68562.6(2.4) keV we find excellent agreement not only 
with AME03 (see Tab. HI) but also with the Penning-trap 
mass measurements at fSOLDE [si and JYFLTRAP 



Rb 



Our measured mass excess for ^''Rb (ri/2 =169.9 ms) 
of -58519.2(2.1) keV lies within the uncertainty of the 
JYFLTRAP value [13, but it differs by 163 keV (5.3cr) 
from the value adopted by AME03 ^2^. We confirm the 
JYFLTRAP value with an improved precision. The pre- 
vious measurements considered in AME03 were based 
on 13 end-point energies from ^^Rb(/3~)^''Sr. The ad- 
justed Q- value, including the TITAN input, is now 
10063(4) keV, compared to the other measurements: 
10020(50) keV [51|] (0.9(T deviationVl0450(30) keV [52| 
(12.9(7 deviation), 10440(60) keV ^ {6.3a deviation), 
and 10462(40) keV [5^ (lO.Ocr deviation). The TP 
TAN mass value validates the first measurement [5l| and 
greatly improves the precision of the Q- value and mass. 



C. ®*Rb (ground and isomeric state) 

This work presents the first Penning-trap mass mea- 
surement of ^^Rb {Ti/2 =114(5) ms). Our measured 
mass excess is -54318.4(3.4) keV and differs from the 
adopted value in AME03 ^ of -54221.6(50.2) keV by 
97 keV or 1.9cr deviation. The previous mass excess was 
determined from the end-point energy of the /3 spectrum 
of 98Rb(/3-)98Sr H, [H with 80.4% weight and from 



the triplet measurement of (^^Rb, ^*Rb, ^^Rb) '55'] with 
19.6% weight. 

A low-lying isomeric state in ^^Rb is proposed at 
286(128) keV [56^ with a half-hfe of 96(3) ms [57]. The 
adopted value in AME03 takes into account not only Q- 
values from /3 end-point energies as mentioned for ^^Rb, 
but also the mass-triplet measurement of reference 15511 
for the ground state combined with ^*Rb™(/3~)^*Sr [53| . 
In an effort to confirm the energy of the isomer, we 
scanned a range of 630 keV (see Fig.[S^), in which we ex- 
pected to observe the isomer. If the isomer were present, 
two dips would have been visible in the resonance curve. 
Only one dip, consistent with the ground state, was ob- 
served; however, the absence of a second resonance does 
not exclude the possibility of an isomer. The strength of 
the resonances depends on the ratio of the population of 
the isomeric and ground states. Yield measurements at 
ISAC j27] indicated the yield of the ground state to be 
20 times larger than the isomeric state. If this ratio were 
observed with MPET, the signal is expected to be the 
black, solid curve in Fig. [HI from which the isomer can- 
not be detected. If, however, there were only five times 
more ions in the ground state than in the isomeric state, 
the isomeric state would be detectable (red, solid curve). 
To guide the eye, a resonance is drawn as if each state 
alone were trapped and measured with a dashed curve 
(gray for the isomer and blue for the ground state). To 
further confirm the energy of the isomer a higher charge 
state, preferably an isoelectronic series of Ar or Ne cor- 
responding to q — 19+ and q — 27+ respectively, could 
be used for the mass measurement to strongly enhance 
the resolving power. This would allow the implementa- 
tion of dipole cleaning of the ground state to enhance the 
resonance of the isomer. 



D. 



94,97,98,99 



Sr 



For all the measured Sr masses the uncertainties were 
reduced. The mass excess of ^*Sr is known to 7 keV accu- 
racy from measurements at ISOLDE [2^ . Our measure- 
ment agrees and improves the accuracy by a factor of 4 
to 1.7 keV as presented in Tab.|T]and Fig.[71 For ^''Sr our 
measured mass differs by 207 keV and 10.8cr deviation to 
AME03 [l^l , but it confirms the mass measurement from 
JYFLTRAP [ll|. In AME03 ^'^Sr is mainly determined 
from ^^Sr(/3~)^''Y. The adjusted Q-value, including the 
TITAN input, is now 7545(8) keV. The previously mea- 
sured Q-value was underestimated: 7452(40) keV [s^ 
(2.3(7 deviation)j 7420(80) keV [53] (1.6cr deviation), and 
7480(18) keV [Hi] (3.6ct deviation). The mass from TP 
TAN is significantly more precise than previous Penning- 
trap mass measurements by a factor of 3. 

The scenario for ^*Sr is similar where our value 
with -66424.5(4.0) keV is in agreement within the er- 
ror of JYFLTRAP ^1 and more precise, but dis- 
agrees with AME03 by 221 keV (8.4(7 deviation). 
AME03 adopted its value from /3 end-point energy ex- 
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FIG. 5. (color online) The cyclotron resonance for «'*Rb^5+ was taken with 80 ms charge-breeding time, 20 ms dipole cleaning 
to eliminate ®*Sr^^+, and 77 ms excitation time in MPET. The solid line is a fit of the theoretical line shape [s^ to the data. 
The left graph shows the resonance off center to include the range of the proposed isomer ^*Rb™ in frequency space indicated 
by a dashed blue square. The graph on the right shows another measurement with a centered resonance. See Tab.Ufor results. 




-40 -30 -20 -10 10 

- 8 700 027 (Hz) 



FIG. 6. (color online) This calculation displays the theoretical 
line shape of the isomer 98|:j^|jmi5+ dashed gray and the 
ground state ^^Rb^*"*" in dashed blue if only one of them is 
present. If both are present the solid lines indicate the line 
shapes for the resonance. If the yield of the isomer is a factor 
of 20 (5) less than for the ground state, the resonance line 
shape is displayed with a solid black (red) line. 



perimcnts, with 95.5% from ^^Sr{p-)^^Y [13] and 4.5% 
from 98Rb(^-)98Sr 

In the case of ^^Sr where the mass potentially plays an 
important role for the r-process and further mass extrap- 
olations, we obtain a mass excess of -62506.8(4.1) keV. 
This value agrees within 2.5a with the JYFLTRAP mea- 
surement [l^, but it disagrees by 321 keV (4. Oct) with 
the mass evaluation ^2^. This is the only case where 
we differ from the JYFLTRAP measurements. Our mass 
excess gravitates slightly away from their measurement 
towards AME03 [2^] where the adopted mass stems from 



/3 end-points from ^^St{P-)^^Y [H, with 91% weight and 
99Rb(/3^)99Sr with 9% weight. 



V. DISCUSSION OF THE RESULTS 



The masses presented in this work impact nuclear 
structure and nuclear astrophysics studies. The atomic 
mass evaluation adopted its final mass value from a 
weighting of mass values obtained from different experi- 
mental techniques. Therefore, links between nuclei exist, 
and our mass values influence other mass values as well. 
This is shown in a more systematic way in Tab.HIl Listed 
here are the most important contributing data that are 
used to determine the mass of the nuclide as it will appear 
in the next mass evaluation 56]. This upcoming mass 
evaluation follows AME03 [221]. It takes into account 
correlations between mass values due to relative mass 
measurements. To determine the evaluated mass value, 
one uses a least-squares method weighted according to 
the precision with which each piece of data is known. 
The measured Rb and Sr mass values not only influence 
their own isotope's evaluated value, but they also influ- 
ence other mass values in the neutron-rich region and 
improve their precision due to links. In various isotopes 
e.g. 9^Zr, ^''Zr, ^"^^Nb, and ^'^'^Nb the mass excess changes 
by one standard deviation. A global mass evaluation was 
carried out and the impact was investigated. However, a 
detailed investigation of nuclei not measured in this work 
lies outside of the scope of this pub lication and can be 
seen in the new mass evaluation 15611. 
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FIG. 7. (color online) Comparison of the mass excesses determined in the present work (red squares), with JYFLTRAP [la. [l7l] 
(turquoise triangles) and in AME03 [l^l (blue circles) for (a) Rb isotopes and (b) Sr isotopes. The inset displays an enlarged 
view. 



TABLE II. Most important contributing data to and their influences on its mass as it will appear in the next mass evaluation [56l | 
following AME03 [22] and extended by the TITAN masses from this work. Influences are given as a percentage. Experimental 
techniques displayed such as (^''Rb vs. *^Rb) indicate Penning-trap mass spectrometry and (^^Rb(/3~)®^Sr) /3 end-point energy 
experiments. 



Nucleus 


Influence 


TITAN 


Influence 


Others experiments 


Evaluated ME (kcV) 


«*Rb 


70.2% 


s^Rb vs. «^Rb 


29.6% 


s^Rb vs. «**Rb [17] 


-68562.3(2.0) 


^^Rb 


87.0% 


^^Rb vs. **^Rb 


12.9% 


»^Rb vs. 8«Rb [17] 


-58518.5(1.9) 


s»Rb 


100% 


3»Rb vs. «^Rb 






-54317.7(3.4) 


«*Sr 


98.4% 


s^Sr vs. «5Rb 


1.6% 


»^Sr(/3-)3*Y [51] 


-78845.1(1.7) 


^^Sr 


87.4% 


^^Sr vs. «5Rb 


12.6% 


^'Sr vs. '''Zx [16] 


-68581.8(3.4) 


8»Sr 


85.2% 


^^'Sr vs. *^Rb 


14.8% 


3»Sr vs. 9^Zr [16] 


-66425.6(3.7) 


^^Sr 


75.9% 


•■'^Sr vs. «5Rb 


24.0% 


s^Sr vs. 99Zr [16] 


-62511.7(3.6) 



A. Nuclear structure findings 

To assess the impact on nuclear structure, Fig. [5] (a) 
shows the mass surface defined by the isotopic two- 
neutron separation energies 82^- It illustrates the de- 
formation for Z = 36 — 40 (Kr to Zr) with data from 
AME03 in open black circles, recent mass spectrometry 
results from JYFLTRAP (for Rb, Sr, Y, Zr, Nb, and Mo 
masses [iMl) and ISOLTRAP (for Kr masses [i,!!^]) 
in filled black circles, and the neutron-rich Rb, Sr data 
from TITAN (this work) in red squares. Extrapolated 
masses are disregarded. 

The smooth trend of the S2n as seen for Kr (see 
Fig. [8] (a, b) and ^) is interrupted for other isotopic 
chains indicating a sudden change in deformation. The 
new data from this work agree with previous experiments 
showing an onset of large deformation for A « 100 nu- 
clei with N > 60. This can be strongly seen for Rb and 
Sr. Previous work indicates a rapidly changing behav- 
ior in nuclear structure in the region of 58 < A'^ < 61. 
To visualize the deformation, the difference between the 



S2n of the isotoncs A^ = 61, A = 59 (black circles) and 
A^ — 63, A^ = 61 (blue triangles) versus the proton num- 
ber is shown in Fig. [8| (b). This observable illustrates 
the so-called quantum nuclear shape transition [6l|, l62j 
between Af = 61toA^ = 59. A shape transition is clearly 
visible for Rb to Mo, whereas Kr presents the lower limit 
and Tc, Ru the upper limits, respectively. With the in- 
put from this work indicated in red squares, we obtain 
new data points extending to more neutron-rich isotopes. 
In contrast, the difference in S2„ for A^ — 63, A^ = 61 
displays the smooth behavior again. The slope of 82™ 
forms its smooth trend again, which clearly indicates a 
strengthening of one nuclear shape. For the Sr isotones 
from A^ = 61 to TV = 59, previous data (AME03, dashed 
black lines in Fig. [8| (b)) showed no signature of unusual 
behavior, while our data and [l6| strongly display the 
shape transition in Sr. 
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FIG. 8. (color online) (a): Two-neutron separation energies (S2n) for Z = 36 — 40 (Kr to Zr) versus neutron number. In open 
black circles data from AME03 [221 is presented. For comparison in filled black circles results from recent mass spectrometry 
experiments from JYFLTRAP and ISOLTRAP [1, H^l are shown extended by the TITAN masses. The TITAN mass 

input itself is indicated by red squares, (b): Differences between the S2n of the isotones A'' = 61, = 59 in black circles, 
N = 63, A = 61 in blue triangles, and previous AME03 [l^l data in dashed lines versus the proton number. This observable 
illustrates the quantum nuclear shape transition between A^ = 61 to A^ = 59 with Kr and Tc, Ru being the lower and upper 
limits. In contrast the difference in S2n for N = 63, A' = 61 displays the normal behavior again. 



B. Astrophysical implications 



The masses measured in this work are relevant for a 
range of different types of r-process models. As an ex- 
ample we explore the astrophysical implications of our 
results with a parameterized, fully dynamic r-process 
model following Freiburghaus et al. The model is 

inspired by the conditions that might be encountered in 
high entropy winds emerging from the nascent neutron 
star in a core collapse supernova explosion. As a start- 
ing point one assumes a fluid element that is heated to 
a very high temperature (T «9 GK) where the composi- 
tion is essentially protons and neutrons, with the electron 
abundance Yg being set by weak interactions. The fluid 
element then undergoes a rapid expansion at constant 
velocity v, Ye, and entropy S {S displayed in entropy 
per baryon in multiples of the Boltzmann constant). We 
choose similar model parameters to Hosmer et al. [ll|, 
i.e. Ye — 0.45 and a velocity v = 7500 km/s. The model 
is coupled to a full reaction network with 5410 isotopes 
that includes all relevant charged particle, /? decay, and 
neutron capture rates that ensue. For unknown masses, 
mass extrapolations from [22*1 and calculated values from 
the finite-range droplet mass model (FRDM) [63r] were 
used. Masses enter exponentially in the calculation of 
(7,n) photo disintegration rates from the forward (n,7) 
rates via detailed balance. Calculations are carried out 
for a grid of entropies with the resulting isotopic abun- 
dances being added up with equal weight. Low entropies 
lead to low neutron-to-seed ratios and a weak r-process 
producing mainly lighter r-process isotopes, while higher 
entropies lead to more extended reaction paths all the 
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FIG. 9. (color online) Calculated r-process abundances 
as a function of mass number summing all entropies for 
AME03high (dashed red) and AME031ow (soUd black) neu- 
tron separation energies for ^^~^^Rb and ®'^~^°°Sr. Also 
shown for comparison are the solar r-process residuals (filled 
blue circles) 64]. 



way to the heaviest elements. An entropy range from 40- 
260 is sufficient to capture all entropies that contribute 
to the r-process. The attractive feature of this model is 
that it is inspired by the conditions one might encounter 
in high entropy winds from nascent neutron stars in core 
collapse supernovae, and that the solar system r-process 
abundance pattern can be reproduced reasonably well 
with just two free parameters - Y^ and v [loj . 

To explore the relevance of the masses measured in this 
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FIG. 10. (color online) Calculated r-process abundances as a function of mass number for a single S = 100 component for 
several data sets of neutron separation energies for ^^~^^Rb and ^^"^'^''Sr. The left figure (a) shows results for variations in the 
neutron separation energies using AMEOShigh (dashed red) and AMEOSlow (solid black); The right figure (b) shows results 
using neutron separation energies from this work varied by the new experimental uncertainties with Scr (NEWhigh and NEWlow 
respectively). In addition in the lower panels, both figures display the relative difference between calculated abundances using 
the high and the low neutron separation energies with |high — low |/ [(high + low)/2]. 



work for r-process simulations we performed two calcula- 
tions where we varied the neutron separation energies of 
96-99Rb and gj^j^gj. (AMEOShigh) or all 

down (AMEOSlow) by their AMEOS [13 3a errors. This 
results in variations of ±40, ±60, and ±140 keV for the 
Rb isotopes, respectively, and of ±SS, ±SS, ±84, and 
±150 keV for the Sr isotopes, respectively. We choose 
a Scr variation as such deviations are not uncommon for 
non-Penning trap mass measurements. In fact, the aver- 
age deviation to the new masses determined in this work 
is Scr, but extends to values as high as 6cr for the neutron 
separation energies of ^^Sr and ^^Sr. 

Fig. |9] shows the resulting composition produced by 
the r-process for both cases. Although the masses used 
for this figure stem from AMEOS 22] with AMEhigh and 
AMElow for the isotopes measured in this work, it clearly 
illustrates that the new masses introduce significant vari- 
ations in the composition around A = 95 — 100 due to 
their deviation from AMEOS of up to 6a in neutron sep- 
aration energies. The affected entropy components are 
about 5 = 70 — 110. The component that is most dra- 
matically affected by the new masses is the S = 100 
component, which is shown in Figs, [l] [TOl a) and (b). 
The reaction flows up to A ss 90 are characterized by a 
complex network of charged particle and neutron induced 
reactions and their inverse. This charged particle process 
provides the seeds for the r-process which then occurs at 
a somewhat later stage when temperatures have dropped 
and charged particle reactions have stopped to operate. 
The remaining free neutrons are then rapidly captured, 
driving the composition to more neutron-rich species and, 
via (3 decays, up to heavier elements. The Sr isotopes are 
located in the transition region between these two types 
of reaction sequences and mark the lightest element in- 



volved in a "rapid neutron capture" reaction sequence at 
this entropy (at higher entropies neutron capture starts 
at lower element numbers, but the r-process converts all 
nuclei into heavier species so there is no longer a con- 
tribution to the A 100 mass region). Higher neutron 
separation energies for the Sr isotopes shift the reaction 
flow towards more neutron-rich nuclei. For low 5„ the 
dominant Sr waiting point, the point in an isotopic chain 
where the /3 decay into the next isotopic chain occurs, is 
^®Sr. For high 5„, a signiflcant fraction of the reaction 
flow proceeds via neutron capture on ^^Sr to the ^"^Sr 
waiting point, leading to an increase in the production of 
A = 100 nuclei. 

With the new Penning trap masses measured in this 
work, combined with the work of Hager et al. [l^, the 
contribution of mass uncertainties in neutron-rich Rb and 
Sr isotopes to the r-process abundance pattern becomes 
negligible. This is shown in Fig. [TUT b) which shows es- 
sentially no change when varying the new neutron sep- 
aration energies within their new Scr uncertainties, with 
NEWhigh being the variation up and NEWlow the vari- 
ation down. 

Also shown in Fig.lHlare the solar r-process abundances 
[G^ I , revealing the common problem of all models of this 
type in reproducing the solar composition of the light 
r-process elements with A < 115. While it is apparent 
that the mass uncertainties of the nuclei considered here 
cannot explain this discrepancy, our measurements are a 
good step towards removing the nuclear physics uncer- 
tainties so one can better characterize the disagreement. 

A somewhat surprising result is that changes in the 
Rb and Sr masses in the A = 9Q — 100 range result in 
significant abundance changes across the entire 5* = 100 
component with significant changes occurring for mass 
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numbers as low a,s A = 70. Clearly the neutron capture 
reaction flow in the Sr region feeds back into the nucle- 
osynthesis of the charged particle reaction sequence. The 
only possible explanation is that the switch in Sr waiting 
points affects the free neutron abundance, which indeed 
is the case here. The shift of the reaction path towards 
more neutron-rich nuclei for AMEOShigh leads to a reduc- 
tion in the neutron abundance. Several conclusions can 
be drawn from this. Firstly, masses can affect the final 
r-process abundances globally, including the production 
of nuclei with lower mass number. Secondly, neutrons 
clearly play an important role in shaping the composi- 
tion produced by the charged particle reaction sequence 
leading to an interplay between r-process and seed pro- 
duction. 



VI. CONCLUSION AND FUTURE 
IMPROVEMENTS 

We have measured the masses of 94,97,98j^|-| ^^^^ 
94,97-99 gj. precision better than 4 keV using the 

Penning-trap mass spectrometer TITAN. This presents 
mass measurements of highly charged ions in the charge 
state q — 15+ and the first direct mass measurement of 
^^Rb. This work also infiuences the adopted mass values 
of several neighboring neutron- rich isotopes [56j . Nuclear 
structure properties, such as the neutron separation en- 
ergy S2n, reveal and validate theory predictions of a sud- 
den onset of large deformation from slightly deformed 
oblate or prolate shapes to strongly prolate shapes in the 
58 < < 61 region for Rb and Sr isotopes. This is 
manifested in Rb and more strongly in Sr with data pre- 
sented in this work. In contrast, the more neutron-rich 



61 < A'' < 63 region reveals no shape transition, and 
the smooth S2n trend is stabilized again. The precise 
and accurate mass measurements of this work will allow 
one to validate theoretical models and refine calculations 
towards more neutron-rich Rb and Sr. This may enable 
models such as the self-consistent mean-field approxima- 
tion based on the DIS- or DIM-Gogny energy density 
functional ([15] and references therein) to be tested un- 
der extreme conditions. 

The differences to previous work of up to llcr deviation 
in mass and 6cr deviation in neutron separation energies 
motivated a study of the r-process. The r-process model 
calculations indicate that the measured masses are now 
known to a precision where their uncertainty does not 
contribute to model uncertainties anymore. 

To strengthen the mass measurement program of HCI 
with TITAN, a preparation trap is presently being built. 
The goal is to increase the precision of mass measure- 
ments with HCI and to compensate for negative effects 
due to the charge-breeding process (efficiency losses, en- 
ergy spread, etc.), and therefore a cooler Penning trap 
(CPET) will be commissioned j65l - l67| in the near future. 
Further extension of mass measurements in this region is 
planned. 
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